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Сравнение методов анализа 
низкочастотных сигналов 


В статье делается попытка сравнения метода микроанализа (МА) с широко применяемым методом анализа 
на основе БПФ. Для адекватности сравнения оба метода разделены на следующие составляющие: объект 
исследования, инструмент исследования, метод исследования. Под объектом исследования понимаются 
совокупности сложных сигналов, полученных с массивов сенсоров. Под инструментом исследования 
понимается набор операций, производимых с объектом, позволяющий разложить сложные сигналы на 
множества элементарных составляющих. Применение метода МА - робототехника и смежные с ней 
области техники. 


Введение 


В работах [1-4] был предложен новый метод обработки сложных сигналов низ- 
кочастотного диапазона — метод микроанализа. Слово «микроанализ» употреблялось 
в контексте того, что объектом анализа мог быть сигнал части наблюдаемого 
процесса и базовые функции (на которые производится разложение) тоже могут 
быть частями сложных функций. Целью введения нового метода было решение 
задач сравнительного анализа сигналов с массивов сенсоров, работающих в области 
низких частот. Предлагалось производить разложение входных сигналов на состав- 
ляющие гармоник - т.н. полупериоды, для определения которых не требовалось вы- 
полнения условия непрерывности на всем интервале существования входного сиг- 
нала. Для каждого полупериода определялись несколько параметров — длительность, 
размах по амплитуде, время существования в общем сигнале, вид базовой гармоники 
полупериода. Для более подробного рассмотрения работы метода по разложению на 
составляющие гармоники были определены основные операции над полупериодами 
и их множествами. 

Отдельно следует рассматривать возможность работы каждого из методов в авто- 
матическом режиме (то есть без вмешательства человека-оператора или исследова- 
теля). Реализация этой возможности открывает возможность работы метода на борту 
мобильного робота. С этой точки зрения вейвлет-анализ слабо применим и не рас- 
сматривается как возможный конкурент. Главная проблема вейвлет-анализа (с точки 
зрения возможности автоматизации процесса) состоит в выборе инструмента, то есть 
в выборе базового вейвлета, с помощью которого обрабатывается входной сигнал. 
На настоящий момент не существует универсальных вейвлетов, дающих одинаково 
хорошие результаты на всех видах сигналов. Не существует общепринятых крите- 
риев выбора вейвлетов — это целиком прерогатива человека-исследователя. Кроме 
того, многие результаты вейвлет-анализа (частотный спектр, фазовый спектр) пред- 
ставляют собой сложную структуру данных, пока плохо поддающихся автоматизи- 
рованному сравнению. 
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| Объект исследования — сигналы с массива сенсоров 


Принимается, что непрерывный аналоговый сигнал с каждого сенсора, входя- 
щего в массив, преобразуется аналого-цифровым преобразователем (АЦП) в отдель- 
ный цифровой поток данных. Метод преобразования — импульсно-кодовая модуляция 
(ИКМ или РСМ в общепринятой аббревиатуре). Шаг времени квантования является 
постоянным и, по крайней мере, на один порядок меньше времени существования 
периода наиболее высокочастотного сигнала. Это означает, что для сигнала с час- 
тотой 4 кГц (период сигнала равен 250 мкС) частота квантования должна быть не 
менее 40 кГц (квант времени равен 25 мкС). Синхронизация отдельных потоков дан- 
ных производится исключительно по времени и ошибка синхронизации не превы- 
шает кванта времени между двумя отсчетами АЦП. Оцифровка данных может произво- 
диться АЦП со встроенным аналоговым коммутатором, но при этом требования к 
шагу квантования по каждому из сигналов остаются прежними. 

Точность работы АЦП — ошибка не более одного младшего значащего разряда 
при одном и том же источнике опорного напряжения для всех АЦП. Этим обеспе- 
чивается минимальный разброс сигналов после работы АЦП. 

Цифровые потоки записываются на носитель в манере «стереосигналов», т.е. 
сенсоры объединены в пары по какому-то принципу и в каждом файле, отождест- 
вленном с двумя потоками сигналов, четные позиции занимают данные с одного 
сенсора, а нечетные -— с другого. Никакого другого смысла в понятие «стереосигнал» 
не вкладывается — просто использован подходящий стандартный формат записи 
данных, для которого имеется множество готовых программ обработки, в том числе 
позволяющих извлекать сигналы отдельного сенсора (канала), образовывать новые 
«стереопары» сигналов. 

Ссылки на отдельные участки сигналов и их цифровые отображения делаются 
либо во временном масштабе (микросекунды, или мкС), либо просто по номерам 
квантов. 


2 Инструменты исследования сигналов в БИФ 


В методе анализа с применением БПФ [1] можно выделить несколько инстру- 
ментов. Это собственно прямое быстрое дискретное преобразование Фурье, с по- 
мощью которого конечная входная последовательность хо, хр, Х2,..,хм! (в общем 
случае комплексных) формулой вида 


Ее ы (1) 


ИЛИ 
МХ == Ве, +Ить 
преобразуется во множество коэффициентов гармоник Х = {Ху,...,Хм и. 
Обратное дискретное преобразование Фурье (ДПФ) заключается в восста- 


новлении последовательности ху, х/, Х»...Хм.1 Из имеющейся последовательности Ху, 
Хь ХЪ,....Хм.1 по формуле: 


1 М-— 72 пп 
=> Же о (2) 
К=0 
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Амплитуда А», частота и и период Т»ь каждой из гармоник связаны с коэф- 
фициентами Х» формулами: 


Ак= ИМх(Ве?, + 1т?у!?; (3) 
и = % (4) 
Ть = ТА. (5) 
Фаза ‹ф, каждой гармоники определяется формулой: 
ф‚ =ага® — : (6) 


Можно говорить, что (1) - (6) — это инструменты БИФ, позволяющие получить 
из входного множества амплитуд сигнала множество гармоник сигнала, обладающих 
своими характеристиками. 


3 Инструменты исследования сигналов в МА 


В методе МА можно выделить несколько инструментов. Это собственно разбие- 
ние входной последовательности амплитуд реального сигнала ху = {Ху Хи, ...Хь.„Хий 
на восходящие и нисходящие реальные полупериоды по условиям: 


х1(1-1) ХИ) <х+1), 0 <1< п, (1) 


х- > ха > ха+1, 0 <1< п. (8) 

При этом каждый реальный полупериод х, приобретает характеристики вида 
образующей функции — //(не определена на первом шаге), размаха амплитуды — ААх, 
длительности по шкале времени -— ат,. 

Метод МА предусматривает наличие базы элементарных полупериодов запи- 
санных в идее матриц вида — М;/, АА, ат, где / — вид образующей элементарной 
функции, АА; — размах амплитуды элементарного полупериода, т; — длительности 
по шкале времени элементарного полупериода. 

Поиском в базе элементарных полупериодов, каждому реальному полупериоду 
ху ставится в соответствие элементарный полупериод 


УЕ МУХ, АА, ату. (9) 


Операцией вычитания амплитуд элементарного полупериода из амплитуд 
реального полупериода получается множество амплитуд разностного сигнала: 


2 = Ху Уь (10) 


Имеется численная оценка качества преобразования в виде суммы разностей 
амплитуд реального и элементарного полупериодов: 


=, В 


В результате операций (7)- (11) из множества амплитуд входного сигнала 
извлекается множество составляющих элементарных полупериодов, каждый из 
которых имеет свою амплитуду, частоту и период (они заданы в матрице М;/, АА, 
ат). Фаза каждого полупериода определяется временем начала реального полупе- 
риода по шкале времени. 

Операции (7) — (11) можно определить как один уровень декомпозиции метода МА. 
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Полученное множество амплитуд 7, можно опять подать на «вход метода» и 
получить на «выходе метода» новые множества у, 2 и К. Процесс декомпозиции 
можно продолжить до тех пор, пока размах амплитуд реальных полупериодов не 
станут меньше допустимого, наперед заданного значения или пока длительности реаль- 
ных полупериодов не станут меньше допустимого, наперед заданного значения. 

По завершению работы метода МА будет получено множество элементарных 
полупериодов У и, поскольку каждый у У обладает характеристикой длительности, 
можно говорить, что У — это множество гармоник сигнала, обладающих своими 
характеристиками. 


4 Сравнение инструментов исследования 


Сравнение инструментов метода БИФ и метода МА можно провести по следую- 
щим параметрам: ограничения на входной сигнал, выбор границ применимости мето- 
дов, вычислительная сложность методов, представления выходных данных. 

Ограничения на входной сигнал влияют на работоспособность методов в целом 
и возможность получения достоверных выходных данных в частности. Метод БИФ 
принципиально не имеет никаких ограничений ни по амплитуде, ни по частоте 
входного сигнала. Метод МА имеет ограничения по амплитуде и частоте входного 
сигнала — это связано с представлением полупериодов. 

Выбор границ применимости методов связан с видом сигналов (гармонические, 
прямоугольные и т.п.). Метод БИФ требует определенной дополнительной обработ- 
ки для негармонических сигналов и дает неверные результаты на смешанных сиг- 
налах. Метод МА целиком зависит от базы элементарных полупериодов, которая 
может пополняться во время его работы. 

Вычислительная сложность метода БПФ обычно определяется как сложность 
различных реализаций алгоритма БДПФ и колеблется от №М02>М до №, где М — число 
обрабатываемых точек. Яркий тому пример — прямая линия. Метод Фурье-анализа все 
равно произведет как минимум М№102>М вычислений, а метод МА обнаружит прямую 
линию на первом шаге и завершит вычисления. Еще пример — сигнал одной частоты — 
Фурье-анализ произведет все №02›М вычислений, метод МА выполнит один уровень 
декомпозиции над всеми полупериодами в исходном сигнале. Кроме того, метод 
БИФ принципиально требует вычислений в формате с плавающей точкой, а для 
реализации метода МА достаточно операций с целыми числами. При одном и том же 
количестве операций (на процессоре не имеющего команд с плавающей арифмети- 
кой) метод МА будет работать примерно в 10 раз быстрее метода БПФ. 

Представление выходных данных важно для хранения результатов и их после- 
дующей обработки. Метод БПФ на любом наборе данных будет всегда иметь фиксиро- 
ванное число коэффициентов, зависящее только от выбора сетки частот, по которым 
идет разложение сигнала. Метод МА принципиально дает различное число выход- 
ных данных, зависящее от сложности и длительности входного сигнала. На простом, 
моночастотном сигнале, метод МА будет давать меньшее число данных, чем метод 
БПФ, а на сложном, длительном сигнале метод МА будет давать на выходе больше 
данных, чем метод БПФ. При этом выходные данные обоих методов — это гармоники 
сигналов. Но для метода БИФ — это гармоники только периодических сигналов, а 
для метода МА - это гармоники любых сигналов, заложенных в базы элементарных 
полупериодов. С этой точки зрения метод МА является более универсальным, т.к. 
позволяет работать с практически любыми видами сигналов и его выходные данные 
могут быть представлены в терминах этих сигналов. 
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Важно заметить, что выходные данные метода БИФ — это гармоники непре- 
рывных функций, а в методе МА - это гармоники любых функций, в совокупности 
представляющие непрерывные функции. 


5 Методы исследования в БФ 


Методы исследования данных, полученных с помощью инструмента исследо- 
вания, основаны на группировке полученных данных и рассмотрении характеристик 
изменения данных. Общепринятыми являются такие характеристики, как Фурье- 
спектр мощности, амплитудный Фурье-спектр, фазовый Фурье-спектр. 

Фурье-спектром мощности называется последовательность: 


=, еб. (12) 


где Р»„— мощность А-й спектральной составляющей. 
Амплитудным Фурье-спектром называется последовательность: 


РЕВ Е (13) 


Фазовый Фурье-спектр определяется из следующих выражений: 


мою 
Ве, (14) 


ф‚ =0,К=0. 


ф‚ = ага 


Как видно, это интегральные характеристики, относящиеся к «окну» данных, 
на которых была построена данная гармоника. Ввиду единственности образующей 
функции, амплитудный и фазовый спектры для каждого сенсора в массиве — это 
строки данных и, соответственно, для массива сенсоров — это квадратные матрицы. 


6 Методы исследования в МА 


В [2], [4] для полупериодов длительностью т; были определены понятия пе- 
риода колебания как Т; = 2хат,, частоты колебания как у; = 1/(2 хат; ) и амплитуды 
колебания как А; = АА(х)). 

В соответствии с этими понятиями, для метода МА можно также можно опре- 
делить спектр мощности для каждой гармоники: 


‚К =0,..К,, (15) 


РЕ = | 


где А- индекс множества полупериодов одинаковой длительности; 
т» — мощность множества полупериодов данной длительности. 

Выражение (15) есть не что иное, как группировка полупериодов одной 
длительности (или одной частоты), которую можно разместить на шкале частот. 
Множество группировок можно отобразить в виде графика амплитудного спектра 
сигнала. Поскольку метод МА имеет множество базовых функций, например, 5т, 
со5, гиперболические, прямая и т.п., то и видов спектров можно получить столько, 
сколько имеется базовых функций. Соответственно можно сравнивать различные 
спектры для одного и того же временного «окна». Таким образом, выражение (15) — 
это множество наборов данных, привязанных ко времени и к базовым функциям. 
Это множество может быть описано в виде квадратной матрицы. 
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Фазовый спектр для метода МА имеет несколько другой смысл — это компо- 
нента ризх; — абсолютное время начала полупериода х; на шкале времени (или индекс 
захвата АЦП). Последнее означает, что метод МА дает дифференциальные харак- 
теристики для каждой гармоники и каждого полупериода данной гармоники. 

Различие базовых функций порождает различие амплитудных и фазовых спек- 
тров и то, что не фиксировалось на спектрах одной базовой функции, может быть за- 
фиксировано на спектре другой базовой функции. Свойство множественности спект- 
ров относится к каждому из сенсоров массива и, соответственно, для одного и того 
же временного «окна» будет получено множество квадратных матриц или, иначе, 
для массива сенсоров результат анализа — это кубическая матрица амплитудных 
спектров и кубическая матрица фазовых спектров. 


7 Сравнение методов исследования 


Метод исследования сигнала на основе БПФ - (12) и метод исследования сиг- 
налов на основе МА - (15) в части амплитудного спектра на одной и той же обра- 
зующей функции дают примерно одинаковый результат, но метод МА не имеет 
«эффекта размазывания», т.к. в нем принципиально не порождаются гармоники, не 
имеющиеся во входном сигнале. 

Точность определения фазы для гармоники в методе БИФ может быть опреде- 
лена только с точностью до «положения окна» + «размер окна». Минимальный 
размер окна для БПФ равна 256 квантов времени получения данных с АЦП. Точ- 
ность определения фазы для гармоники в методе МА это +1 квант времени получе- 
ния данных с АЦП. Таким образом, метод МА дает более точный результат. 

Фазовый спектр в методе БИФ — это результат вычислений и повышение 
разрядности АЦП и уменьшение шага квантования никак не скажутся на этих ре- 
зультатах. Наоборот, фазовый спектр в методе МА - это результат измерений и по- 
вышение разрядности АЦП и уменьшение шага квантования приведут к получению 
более точных результатов (до определенного порога, разумеется). 

Результат анализа для массива сенсоров с применением БИФ -— это квадратная 
матрица амплитудных спектров и фазовых спектров, а с применением МА - это 
кубическая матрица. То есть при одних и тех же исходных данных имеется сущест- 
венно различное количество выходных данных. 


Выводы 


Метод МА имеет ряд преимуществ перед методом БИФ в случае работы с 
произвольными, сложными сигналами. 

Первым преимуществом является большая скорость обработки сигналов. 

Вторым преимуществом является точность определения амплитудного спектра 
сигнала — в методе МА она гораздо выше за счет работы только с гармониками 
(полупериодами), содержащимся в исходном сигнале, а не с заранее определенным 
набором гармоник, как в методе БПФ. 

Третьим преимуществом является точность определения фазового спектра 
сигнала, которая в методе МА гораздо выше за счет привязки гармоник (полупе- 
риодов) исходного сигнала, к абсолютной шкале времени, а не ко временному 
«окну», как в методе БИФ. 

В-четвертых, метод МА позволяет в автоматическом режиме определить вид 
входного сигнала и, соответственно, принять решение о способе его дальнейшей об- 
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работки вплоть до проброса участка входного сигнала. Это совершенно новое свойство, 
которого нет ни в одном другом методе анализа сигналов. 

Метод МА имеет значительно более сложную структуру выходных данных, чем 
метод БПФ. Это можно отнести к недостаткам метода МА. 

Метод МА не имеет строгой доказательной базы именно из-за принципиально 
неограниченного множества элементарных полупериодов. Это конечно недостаток 
метода. Собственно такой же недостаток и у вейвлет-анализа. 
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С.0. Полвцев 

Порвняння методв аналзу низькочастотних сигнал! в 

У статт! робиться спроба пор1вняння методу мкроанал1зу (МА) з широко вживаним методом анализу на 
основ! БИФ. Для адекватност! порвняння обидва методи роздлено на наступн! складов: об’ект дослдження, 
1нструмент дослдження, метод досллдження. Шд об’ектом дослдження розумпюоться сукупност! складних 
сигналив, отриманих з масивйв сенсор!в. Шд 1нструментом досллдження розумпеться наб1р операщй, 
вироблюваних з об’ектом, що дозволяе розкласти складн! сигнали на безл!ч елементарних складових. 
Вживання методу МА - робототехнка 1 сум жн! з нею област! техн ки. 


5.А. Ройл5еу 

Сотрат150п оЁ тео@$ оЁ Ве апа1у$ 5 оЁ10у-Вгедиепсу $101а[5 

ш агасе абетре оЁ сотраг150п оЁа тефо4 оЁ Фе плсгоапа[у$15 (МА) у улдейу аррПе4 тефо4 о{ ®е 
апа[уз1$ оп Фе Баз1$ о ЕЕТ Бесотез. Еог адедиасу оЁ сотраг15оп, Бо тефо4$ аге ату14е4 шо ЮПо\пе 
сотропеп: оБ]ес{ оЁ гезеагсВ, фе гезеагсЬ 1юо01, а гезеагсН плефо4. ТБе об]есё оф гезеагсВ 15 ипдегзюо4 а5 
5е5 оЁ ЧЁИсий 512па15 оЁ Фе зепзог сопго[5 гесеуе4 Нот Шез. ТБе гезеагсЬ {001 15 ипаегзюо4 аз а $её оЁ 
Фе орегайоп$ таде у оБ]есё, аЦо\лпе, №0 зргеаЯ опё ЧР си Е 5121а1$ 10 зе оЁ еетешагу сотропепб. 
АррПсайоп оР тео МА - а гобойсз ап4 а4асепй \уйВ 1 оР1есВитс$ агеа. 


Статья поступила в редакцию 25.02.2010. 
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